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Abstract. We consider 2- and 3-dimensional cubic monocrystalline and polycrystalline 
materials. Expressions for Young's and shear moduli and Poisson's ratio are expressed in 
terms of eigenvalues of the stiffness tensor. Such a form is well suited for studying properties 
of these mechanical characteristics on sides of the stability triangles. For crystalline high- 
symmetry directions lines of vanishing Poisson's ratio are found. These lines demarcate regions 
of the stability triangle into areas of various auxeticity properties. The simplest model of 
polycrystalline 2D and 3D cubic materials is considered. In polycrystalline phases the region of 
complete auxetics is larger than for monocrystalline materials. 



1. Introduction 

In our recent papers we considered mechanical characteristics of 3D crystalline structures of high 
and middle symmetry [lj and of all 2D symmetry systems [2|. In particular we derived explicit 
expressions for Young's E and shear G moduli as well as for Poisson's ratio v depending on 
directional cosines of angles between directions n of the load and the direction m of the lateral 
strain, respectively, and directions of the crystalline symmetry axes. 

Our particular attention was paid to 3D and 2D cubic structures [3], [lj. All initially 
unstrained quadratic and cubic materials are characterized by three parameters si, S2 and S3 
belonging to a half-infinite (s\ > 0) prisms with (stability) triangles (ST) in the base lying in 
the (s2,S3)-plane. 

Among mechanical characteristics, Poisson's ratio v is particularly interesting because the 
phenomenon of a negative v, i.e. solids expanding transversely to applied tensile stress, is 
counterintuitive. Elastic materials with negative v are termed auxetics. Materials with negative 
v for all pairs of vectors n, m are complete auxetics. Materials are auxetics if there exist pairs 
of n, m for which v < 0. Materials are non-auxetics if v > for all pairs of vectors n, m. In 
papers [3], [1] we established regions of ST in which crystalline cubic and quadratic materials 
are completely auxetic, auxetic and non-auxetic. One may expect that complete auxetics are 
particularly promising from the standpoint of technology. 

Many technically important materials are polycrystalline. In view of potential applications of 
crystalline auxetic materials, the question how polycrystallinity influences the auxetic properties 
deserves some attention. 

The elastic properties of polycrystalline materials depend on the single-crystal stiffnesses 
(elastic constants) of the crystallites which build up the polycrystal and on the manner in which 



crystallites are connected. In most cases exact orientations, shapes and connections are not 
known. These gaps in knowledge can be overcome with the help of orientation and grain shape 
distribution functions. Generally, one has to resort to more or less realistic assumptions (cf. 
[5]). We assume that the grains have the same shape and volume and that their orientations 
are completely random. 



2. The stability conditions for initially unstressed crystalline quadratic and cubic 
media 

Mechanical properties of crystalline materials depend on components of the stiffness tensor C 
or on the components of the compliance tensor S. These tensors are mutually inverse, i.e. 

CS = SC = I 4 , (1) 

where UQ^jy) = {&ik&jl + dijdjk) /2. The product of two fourth rank tensors A and B has 

components (AB)^ = Y^, s =iAj,rsBrs,kl- 

In this paper we consider only 2D and 3D cubic crystalline or polycrystalline materials. For 
mechanically stable initially unstrained materials both C and S tensors are positive, therefore 
their eigenvalues cj, sj (I = J,L,M) are positive too ([T], [2], [B]) 

d > 0,s/ > (I = J,M,L). (2) 

Because of Eq. ([1]) sj = cj 1 vanishing eigenvalues cj (I = J,L,M) signal instabilities related to 
phase transitions [7]. 

In the case of quadratic materials [2] 

cj = Cn + C12, cl = 611 — C12, cm = 2C66- (3) 

In the case of cubic materials pQ, [S] 

cj = C n + 2C12, c L = C u - C 12 , c M = 2C 66 . (4) 

In Eqs. ([3]) and ([!]) we used the familiar Voigt's notation 0: Cn = C12 = Cn ; 22 and 

Cm = 612,12- 

For cubic materials a more reasonable choice of independent parameters was proposed by 
Every [9] and for quadratic materials in our paper [1] 

si = 6n + {d- 1)6 6 6, s 2 = , 

si 

6*11 — 612 — 2666 / c \ 

S3 = , (5) 

si 

where d = 2 for 2D, d = 3 for 3D materials. 

We rewrite expressions ([3]) and (JH) using Every's parameters s\, s 2 and S3. In the case of 
quadratic materials we get 

cj = si (2s 2 - s 3 ) , c L = si (1 - s 2 + s 3 ) , s M = 2si (1 - s 2 ) , (6) 

whereas in the case of cubic materials 

cj = y (10s 2 - 6s 3 - 1) , c L = 2^ (1 - s 2 ) , s M = 2|- (3s 3 - 2s 2 + 2) . (7) 



For both cubic and quadratic materials the first of the inequalities ([2]) define stability regions 
in the form of semi-infinite prisms (s\ > 0) with two distinct triangles (STs) at the base in s 2 , 
S3-plane. 

si >0,s 2 < 1, (8) 

2D : (1 - s 2 + a 3 ) > 0, (s 2 - a 3 /2) > 0, 
3D : (1 + 2s 2 ) > |4s 2 - 3s 3 - 1| , (10s 2 - 6s 3 ) > 1. (9) 
Both stability triangles are lying in the (s 2 > S3)-plane (Fig. [1]). The sides of the stability triangles 

a) b) 




Figure 1. a) Division of the stability triangle for quadratic materials into areas of various auxetic 
properties, b) The shaded area is the region of complete p oly cry st aline auxetics. Locations of 
results for graphene [12], [13] and graphene-like ZnO [H] are indicated. 

(STs) are defined by lines 

2D : 4 2D ' X) = 2s 2 , 4 2D < 2) = s 2 - 1, 4 2D ' 3) = 1, 

3D : 4 3A1) = W3 - 1/6, 4 3D ' 2) = W3 - 2/3, sf d>2) = 1. (10) 

A point of the stability triangle together with the value of s\ represents an existing or 
potentially existing crystalline material. For acoustic and elastic characteristics the parameter 
si is a scaling parameter (cf. [7], [TO]). In Figs Q] and [2] we show changes of the division of STs in 
two regions into various auxetic properties. Points of area with lighter shade represent complete 
monocrystaline auxetics, whereas in the area with darker shade lie complete polycrystaline 
auxetics. 

For symmetry systems other than cubic and isotropic, as well as for the oblique symmetry 
system, elastic properties are defined in spaces of dimensions higher than 3. In the case of 
rectangular symmetry system the space of parameters can be effectively reduced to a 3D space 
(cf. Appendix). Monocrystalline stable isotropic materials are located on the interval of the 
lines defined by equation S3 = belonging to the appropriate triangle. 




Figure 2. a) Division of the stability triangle for cubic materials into areas of various auxetic 
properties, b) The shaded area is the region of complete p oly cry st aline auxetics. Locations of 
several classes of elastic materials are indicated (cf. |17j). 



3. Mechanical characteristics of quadratic and cubic materials 

Mechanical properties of an elastic body are characterized by Young's and shear moduli and 
Poisson's ratio. Recent interest in crystalline materials has made the understanding these 
characteristics increasingly important. In papers pQ and [2] we derived expressions for them 
which do not depend on the choice of Cartesian coordinate system. The Young modulus E{n) 
depends on the (direction) cosines that unit vector n of the load makes with the crystalline axes. 
The shear modulus G, as well as the Poisson's ratio u), additionally depend on the direction 
cosines of the vector of lateral load m. 

Expressions for these characteristics can be written in terms of the eigenvalues cj, cl and cm 
For quadratic materials we obtain 

771 / S 2CJC L C M 1 s 

E l( n ) = ; — ^2—, i"£> ( n ) 



clCm + cjc M (nj - nl) 2 + ^cjc L n\n\ 

1 c L c M 

2 cm {mini - m 2 n 2 ) 2 + cl (min 2 + m 2 ni) 



G q (n,m) = ±— ; (12) 



_ c L c M + cjc M {ml - ml) (n| - n|) + 4:Cjc L mim 2 nin 2 , , 

clcm + cjc M {n\ - nl) + ^cjc L n\n\ 



For cubic materials we get 



3cjc L c M 

{clcm ~ cjc L + 3cjc M ) + 3cj {c L - c M ) Yh=i n \ 



^c(n) = — — ; 4 , (14) 



G c (n, m) = — ; — 2-j, (15) 



c.jc L c M 

cjc M + 2cj (ex - c m ) J2t=i mfnj 

, c L (cm - cj) + 3cj (cl - cm) E|=i "^n, 2 
- j/ c (n, m) = — o -. (lb) 

[c l cm ~ cjc L + Scjc M ) + 3cj (cl - Cm) Ei=l n i 

Note that expressions (I11H16P allow one to calculate E, G and v for a given pair of direction 
vectors (m, n) and for chosen values of si, S2 and S3. In the case of polycrystalline materials 
numerical calculation provides one with analogous results. Performing calculations for a mesh 
of points of STs we obtained the appropriate maps [T3]. These maps were calculated for s\ = 1. 

The direction vector of applied tension n and direction vector m in which the lateral 
expansion/contraction is measured are mutually perpendicular, i.e. nm = 0. For quadratic 
materials 



ni = ni2 = cos tp, 

ri2 = —mi = sin 93, (17) 



where < ip < 2tt. 

In the case of cubic materials |11| 



n\ = cos a cos ip cos 6 — sin a sin 9, 
112 = cos a cos ip sin 8 + sin a cos 9, 

713 = — cos a sin (18) 



and 



m\ = sm</?cos#, 
m2 = sin (p sin 9, 

ms = cos tp, (19) 

where 

-7r<a<vr, 0<c/?<7T, 0<9< 2vr. (20) 
4. Mechanical characteristics on the sides of stability triangles 

For directions of high symmetry the mechanical characteristics can vanish, signaling phase 
transitions (cf. [7]). New phases are indicated in Fig. [TJ 

Consider quadratic materials. For n^ = (11), = (11) 

2cjc M r _c L _ cj-c M , , 

&q - ; ) G 9) — TT-i v q- : • \ lV ) 

CJ + C M 2 Cj +C M 

This means that for E q (ni,mi) vanishes on the top and on the vertical side of ST, whereas 
G q (ni, mi) vanishes on the bottom side of it. 
For n 2 = (10), m 2 = (0 1) 

E q = ^,G, = f,v q = C -^±. (22) 
cj + c L 2 cj + c L 



For n 2 and ni2 E q vanishes on the top and bottom side of ST, whereas G q vanishes on the 
vertical side. 



For n 3 9) = (4 9 ' 1} njf' 2) ) and m 3 9) = (n^ n 3 q ' 1] 



Acjc L c M „ c L c M cj (c L + c M ) - 2c L c M /ooA 

h/q — j (jq — -i Vq — 7 : , V^"J 

*C L CM + CjCL + CjC M C L + C M Cj {C L + Cm) + 2c L C M 



_ _^/^ „(ff,2) _ ;^/^ r„, io _„ ang that E q (n^) 

d vertical side of 
,( c ) _ /n n i\ ™( c ) 



where ra 3 = cos(7r/8) and n 3 = sin(-7r/8). This means that E q (n s q ) vanishes on all sides of 
?«(n<*\m?' 



ST, whereas Grg(n 3 9 \ mi ) vanishes on the bottom and vertical side of ST. 



Similar results are valid for cubic materials. For = (0 1), =(110) Eqs. (|14|) . 
(JTSJ and (HBP give 

n (<0 __(c) . p ^ 3c J c M r _ c L cj-cm , , 

n l > m i • h c - ■z . , Lz c - — , v c - - ■ , (24) 

2c L + c M 2 2cj + cm 

hence i£ c ( n i ) = on the top and bottom side of ST, G c (n^ c \m^) = on the vertical side of 
it. For n 2 c) = (110) , mj> c) = (110) we obtain 

„(c) m (c) . p 3c J c ^ c M 

Ilo 111 , . Fj r — ; ; — , 

2 2 cjc l /2 + 3cjcm/2 + c l cm/2' 

C ~ 2' " c ~ cjc l /2 + 3cjcm/2 + c l c m ' 

(c) (c) 

For and Young's modulus vanishes on all sides of ST, whereas the shear modulus 
vanishes on the bottom side of it. Finally, for rig =(111), mj^ = (011) 

2cj + c L 2c L + c M c L - 2cj 

(c) (c) 

For n 3 and m 3 Young's modulus vanishes on the top and vertical side of ST, the shear modulus 
vanishes on the vertical and bottom side. 

For all considered high symmetry directions Young's and shear moduli are non-negative. 



5. Lines of vanishing Poisson's ratio 

For a chosen direction n of load and lateral strain m the general expressions presented in Sect. [3] 
allow us to find the lines of vanishing v that divide the stability triangles into regions of various 
auxeticity properties. These lines are particularly interesting for n and m directed along high 
symmetry directions. 

Consider the quadratic materials. For the following pairs (n, m) the condition v q (m, n) = 
leads to equations 

(n{ q \ m[ q) ) : s 3 (s 2 ) = 2 (s 2 - 1) , 
(n^, m^) : s 3 (s 2 ) = 3s 2 /2 - 1/2. (27) 

For cubic materials we get 

(nf\ m< c) ) : s 3 (s 2 ) = As 2 /3 - 1, 
(4 C) , r4 c) ) : s 3 (s 2 ) = \s 2 - A ± J^68 S | + 20 S2 -7, 

(n 3 c) , m 3 c) ) : s 3 (s 2 ) = 2s 2 - 1. (28) 

All the above lines of vanishing v are shown in Fig. [TJ 

In our papers [3j, H we introduced the division of the stability triangle into regions of 
complete auxeticity, auxeticity and non-auxeticity. The above regions are indicated on diagrams 
in Fig. [H Inspecting them one can distinguish lines of vanishing v studied in this section. 



Mechanical characteristics of 2- and 3-d cubic polycrystalline materials 

Following the assumptions made at the end of Sect. [T] to find the mechanical characteristics 
of polycrystalline materials one should calculate mean values of expressions defining them. 
Consider quadratic materials and a function F of the angle (p, then in agreement with definition 

(F) = ^ J^dtpFfr). (29) 
Calculating (E q ), (G q ) and (v q ) we obtained 

(E,) = 4CJC1C " 2 ■ (30) 

[2c L c M + cj (c L + c M )} - cj (c L - C M ) 

(G q ) = (31) 

(u q ) = 1 - 4CLCM , (32) 



[2c L c M + cj (cl + c M )] - cj (c L - cmY 

Table [U summarizes the properties of mean values of the mechanical characteristics and 
eigenvalues cj (I = J, L, M) on sides of ST. These properties align with calculated maps of 
(E q ), (G g ) and (v q ) [E]. 



Table 1. Values of (E q ), (G q ) and (v q ) on the sides of ST. J - top side, M - bottom side, L - 
vertical side; U = J, L, M. 



Side 


cu 


(E q ) 


(Gq) 




J 








^Jcmcl/2 


-1 


L 











1 


M 











1 



In the case of quadratic polycrystalline materials the condition (v q ) = leads to equations 
of line 

3g| - 10s 2 + 3 

83 (S2) = 3, 2 - 5 ' (33) 

demarcating the stability triangle to the regions of complete auxetics and non-auxetics. In 
case of cubic materials such a line is found numerically. Inspecting Fig. [2] we note that for 
polycrystalline 2D and 3D cubic materials the area of region of complete auxeticity is larger 
than for monocrystalline materials. 

In the case of cubic materials and a function F of angles a, ip and 8 (cf. Eqs. (|18p and (|19p ) 
the mean value (F) of F(a, if, 6) is defined as 

(F) = -^f + da [ dp sin p [ dOF(a,<p,9). (34) 
ovr^ J-t, Jo Jo 

Since we consider the model of randomly oriented crystallites, the averaged mechanical 
characteristics are isotropic. This means, that if (v) > 0, one deals with isotropic non-auxetics, 
whereas if (v) < 0, a polycrystal is a complete isotropic auxetic. We conclude that in the case 
of considered models of polycrystals the division of the stability triangles simplifies. 



6. Concluding remarks 

We have shown that mechanical characteristics of elastic materials expressed in terms of 
eigenvalues of the stiffness tensor are well suited for the study of their general properties. The 
simplest model of polycrystalline 2D and 3D cubic materials indicates that some monocrystalline 
auxetics become complete auxetics in the polycrystalline phase. 

Appendix 

Elastic properties of rectangular materials are characterized by four stiffnesses Cn, C 22 , Ci2> 
Cqq (cf. [2] and [IS]). Initially unstrained rectangular materials are stable if 

C11 > 0, C 22 > 0, C 66 > 0, C n C 22 > C\ 2 . (35) 

Introduce four new parameters 

si = (Cn + C 22 )/2 + C 66 , (36) 

s 2 = [(C11 + C 22 ) /2 - C 66 ] (37) 

S3 = [(C11 + C 22 ) /2 - Ci 2 - 2C 66 ] /s 1; (38) 

s± = {C 11 -C 12 )/2s 1 . (39) 

Because of inequalities ([36]) - ([39]) 

si > 0, -1< s 2 < 1, (40) 

S4 can be positive and negative. 




Figure 3. Rectangular materials: mechanical stability region in (s 2 , S3, S4)-space. 

Consider 3D space of parameters s 2 , S3, S4. In planes s 2 = a (—1 < a < 1) the last of 
inequalities ([35]) takes the form 

sl + (s 3 - di) 2 < 4, (41) 

where d\ (a) = (3a — 1) /2, d 2 (a) = (1 + a) /2. Inequality (fldl) defines a section of the oblique 
cone with the apex located at the point (s 2 = —1, S3 = —2, S4 = 0). For —1 < a < 1 the cross- 
sections of the cone perpendicular to s 2 axis are circles with radii equal to d 2 {a) and centers 
(s 2 = a, S3 = di(a)) , S4 = (Fig. [3|). For S4 = the cone (jlT]) reduces to ST defined by Eq. 
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